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-  amplitude  of  free  waves 

-^coefficients  of  the  free-waves  resulting  from  the 
h  second  order  free-surface  correction 

-  as  above  for  the  body,  correction 

-  huler  constant  .• 

-  drag  coefficient  ... 

-  drag  coefficient  .. 

-  wave  resistance  (I)  =  O'g/pU'4) 

-  first  order  wave  resistance 

-  second  order  wave  resistance 

-  wave  resistance  resulting  from  second  order 
froe-surface  and  body  effects,  respectively 

-  exponential  integral  ... 

'*  complex  potential  (f  '*»  f'g/U’3) 

-  first  and  second  order  complex  potentials,  respectively 
>  acceleration  of  gravity 

-  depth  of  submersion  (h  =  h'g/tl'2)  • 

-  distance  between  the  sources  j  and  k 

-  forebody  or  body  length  (1,  *  1. 

-  second  order  pressure 

-  strength  of  source  j 

-  thickness 

-  thickness  change  at  source  j 
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t  =  t'/h'  -  thickness  (dimensionless  with  respect  to  h 1 ) 

x',yf  -  cartesian  coordinates  (x,y  =  x'g/U'2,  y'g/U'2) 

z'  =  x'  +  iy'  -  complex  variable  (z  =  z'g/U'2) 

u',v'  -  velocity 'components  (u  =  u'/U',  v  =  v ' /U  * ) 

U *  -  unperturbed  velocity 

w'  =  u'  -  iv*  -  complex  velocity  (w  =  w'/ll') 

a  -  real  part  of  w 
3  -  imaginary  part  of  w 

-  auxiliary  function 
<5^  -  auxiliary  function 

e.  =  2uq!g/tJ'3  -  strength  of  source  j  (dimensionless) 

J  J 

e  »  T'/J<*  -  slenderness  parameter 

1# 

(f>*  -  potential  (<{>  =  (fr'g/tJ'3) 

X „  auxiliary  function 

<l>f  -  stream  function  *  i^'g/lJ'3) 

r,  -  complex  variable 

o'  -  free » surface  elevation  (n  ■=  n'g/U*2) 

«(0  «■  e’^EiCic).  -  function  associated  with  waves 
X,p  -  complex  variables 
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ABSTRACT 

The  second  order  wave  resistance  of  an  arbitrary  distribution  of  sources, 
representing  a  thin  symmetrical  two-dimensional  body,  is  derived  in  a  closed 
analytical  form. 

The  result  is  applied  to  a  few  shapes:  an  isolated  source  representing 
semi-infinite  blunt  body,  an  open  body  with  a  leading  edge  of  a  fine  shape  and 
a  closod  body  generated  by  a  source  and  a  sink.  In  each  caso  the  first  order 
vave  resistance  as  well  as  the  body  and  free- surface  second  order  corrections 
arc  derived  separately. 

The  relevance  of  the  results  to  ship  wave  resistance  is  disouns®d  in 
qualitative  terms.  . 
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I.  INTRODUCTION 

The  wave  resistance  of  bodies  moving  beneath  or  at  a  free-surface  is 
generally  computed  by  linearizing  the  equations  of  flow,  i.e.  by  assuming  that 
the  body  is  thin  or  slender.  This  way  the  free-surface  condition  is  linearized 
and  the  body  is  replaced  by  a  distribution  of  singularities  of  known  strength. 

Till  recently  only  the  first  order  wave  resistance  has  been  computed  in 
most  cases.  Nonlinear  effects  associated  with  the  body  condition  have  been 
investigated  by  Havelock  (192b)  and  (tiesing  and  Smith  (1971)  for  two-dimensional 
flows  and  by  Gadd  (1970) and  others  for  three-dimensional  flows  past  ships.  In 
all  these  works  the  free-surface  condition  has  been  kept  in  its  first  order 
version.  The  procedure  is  not  consistent  in  principle  because  body  and  free- 
surface  effects  are  of  the  same  order  of  magnitude,  in  an  asymptotic  sense  at 
1  east . 

A  complete  second  order  numerical  computation  of  the  wovo  resistance  has 
been  carried  out  by  Tuck  (1965)  in  the  case  of  a  submerged  cylinder.  The  free- 
surface  correction  turned  out  to  he  larger  than  the  hotly  correction  in  this 
case.  Hagan  (1971)  has  shown  that  as  the  cylinder  approaches  the  free-surface, 
nonlinear  effects  have  to  be  incorporated  already  in  the  first  order  approxi¬ 
mation.  Sttlvesen  (1969)  has  determined  numerically  the  second  order  wave 
resistance  of  a  submerged  symmetrical  hydrofoil.  The. body  correction  appeared 
to  he  the  largest  in  this  case,  apparently  because  of  the  important  role  played 
by  the  circulation  associated  with  the  Kntta-Juvkovsky  condition  at  the 
trailing  edge. 

The  computation  of  the  nonlinear  wave  resistance  for  two-dimensional 
bodies  is  of  a  limited  interest  in  applications.  Jin  fact,  the  main  purpose  of 
the  works  mentioned  so  far  was  to  draw  conclusions  on  ship  wave  resistance  from 
the  two-dimensional  results,  at  least  on  a  qualitative  basis.  The  two-dimensional 
solution  is  considerably  simpler  than  that  of  three-dimensional  flow  and  it  is 
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worthwhile,  therefore,  to  explore  first  the  two-dimensional  case  in  order  to 
acquire  insight  in  the  problem  and  experience  in  handling  the  mathematical 
difficulties. 

The  present  work,  dealing  with  the  ws'A.  resistance  of  a  two-dimensional 
source  generated  body,  is  a  continuation  of  Tuck  and  Salvesen  works,  but 
differ  from  them  in:  (i)  the  solution  is  carried  out  for  bodies  of  arbitrar) 
thickness  distribution  rather  than  particular  cases,  (ii)  the  second  order 
wave  resistance  is  obtained  in  a  closed  analytical  form,  rather  than  by  tedious 
numerical  computations, such  that  a  general  analysis  of  the  results  is  possible, 
and  (iii)  the  wave  resistance  is  computed  for  two-dimensional  bodies  resembling 
ships,  i.e.  elongated  bodies  of  a  fine  form  and  without  circulation.  The 
cylinder  has  an  extremely  blunt  shape  and  in  the  case  of  the  hydrofoil  the 
contribution  of  the  circulation  has  no  three-dimensional  counterpart. 

II.  EtvUATltXHS  AND  BOUNDARY  CONDITIONS 

Ke  consider  a  discrete  distribution  of  sources  of  strength  qj  (Pig.  ), 
located  at  x*  »  r* .  y*  «.  0  {j«l,2M.,,n).  The  sources  are  assumed  to  be  at 
the  same  elevation  because  any  thin  body  without  circulation  degenerates  into 
such  a  distribution.  It  is  quite  easy,  however,  to  extend  the  approach  pre¬ 
sented  here  to  any  distribution  of  sources  or  other  singularities. 

be  make  the  variables  dimension less  as  follows 

x  «  x*g/irs,  y  *  y'g/ti’^,  ’t*  x*iy-  *  sVg/U,s,  n*  n’g/U‘a 

w  *  u-iv  “  w'/OV,  f  *>  #  ♦  *  f’g/0'3,  h  *  h'g/W**,  (1) 

c.  *  istq'g/D1'5  ' 

J  J 
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whore  f'  is  the  complex  potential,  w'  is  the  complex  velocity,  h'  the 
submergeancc  depth,  U'  the  unperturbed  velocity  and  n'  the  free-surface 
elevation  (Fig.  1). 

If  wo  relate  qj  to  the  change  in  the  thickness  Tj  of  a  body  moving 
in  an  infinite  flow  domain,  then 


Tjft 


(2) 


f  and  w  are  uow  expanded  in  asymptotic  series 


t\  nil 

f  a  s+r  ♦f  ♦... 


(3) 


W  a  £;•  *  1  ♦  W1  ♦  ♦ 


these  being  the  only  terms 


where  f*»  wf  »  0(e)  and  f^,  *»  t 

considered  herein. 

The  analytical  functions  f  (?)  and  f  (s)  satisfy  the  *V)  lowing  . 
linearised  tree* surface  conditions  (see,  for  instance,  hehausen  and  Lai tone, 
1360)  • 


♦  if1) p*(x)  «  (» 

*  s 


(y  *  U) 


V. 


u 


»* 


Ji 


it 


.3  ir.  1*2  *  1 


♦.  if  )  *  !»*•(*)  *  v  S«T  >  s  Cv n  ■  ■  if  *  U) 


(5) 

(6) 


and  the  radiation  com!  it  inn 


fj  '»  ®  , 


(x  ■*  -  «■) 

Use  free-surface  elevation  n  e  *  *  . . .  ^  given  by 
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II,  n 
n  (x) 


-  il»*(x,h) 

n1  Cx)2  ,11,  ,, 
— k ~ —  -  #  (x,h) 


(S) 

(9) 


The  body  condition  reuuces  to  the  requirement  that  at  emu  order  the  thick- 
chan£ 
stipulate 


ness  change  at  the  source  j  i  equal  to  Tj .  At  second  order  we,  therefore. 


e . 


^  =  2  +  2tt  in(z“x-)  *  °(e3) 
^  '1 


Z  -  X. 


j 


•  e . 
J 


(10) 


The  purpose  of  the  following  sections  is  to  determine  th.-  wave  resistance 
of  the  source  distribution. 


III.  THE  FREE  WAVES  BEHIND  A  SOURCE  DISTRIBUTION 

To  compute  the  wave  resistance  we  need  only  the  expression  of  the  free 

waves  trailing  far  behind  the  body.  Moreover,  the  solution  for  a  pair  of 

sources  may  be  easily  extended  to  an  arbitrary  number  of  singularities.  For 

tho  s'ke  of  simplicity  wo  carry  out  first  the5  analysis  for  two  sources,  say 

of  strength  t:  ■  and  e.  ,  located  at  x.  and  x,  Y  respectively . 

J  K  J  K 

1,  First  order  solution 


Tho  first  order  .solution  is  well  known  (Wehausen  and  Laitone,  1960) 
and  is  given  hero  for  convonionce  of  reference 


A 


c.  e.  e. 

fjj,  »  ’j'-  4n(2-x^)  +  27  &n(z-Xj.)  +  ^7  £n(z-Xj -2ih)  + 

+  ~~  £n(z-xk-2ih)  -  «(*.-Xj-2ih)  -  w(z-xk«2ih) 


(U) 
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J  _ _ _  • 

ii  2tt  z-x.  2n  z-x,  2n  z-x.-2ih  2v  z-x,  -2ih.  , 

J  h  j  k 

ie.  ir, 

*  u(z-x^.-2ih)  +  ~r  w(z-xk~2ih)  '  (12) 

The  function  w(0  defined  as 

.  _ •  •  .  t,  i\ 

'  oj(0  =  _a  +  iB  =  e"Uf.i(U)  =  “e"iC  /  ~-dA  ■  ,  (13) 

‘  •  '  -00  \  _v  •  .  *  “  .  •  '•  ■  .  v  ' 

represents  the  trailing  singularities  associated  with  the  free  waves.  The 
properties  of  w(s),  which  will  be  extensively  used  in  the  sequel,  are 
discussed  in  Appendix.  • 

The  far  downstream  potential  is  obtained  by  expanding  (11)  (see  Hq.  (65)) 

T  -i(z-x.-2ih)  :  -i(z-x. -2ih) 

.  :  'f!k -2ie.  e  J  -  2iek  e  k  (x  +  «)  (14) 

and  the  fxvo  waves  profile  (by  liq .  (8))  is  given  by 

ix  *  -  ■  ’  ix 

nj^(X)  “ .  ■I-B4.-2i)e":l(6j_©  *  +  eke  k)e'lx1=  ■  ■"> 

=  2e  ^((e.cbs  x.+e.cos  x,  )cos  x  +  (e.sin  x.+e.sin  x.)sin  x] 

'  •  ‘  J  1  .8  '  A  J  J  ..  A-  K 

-  '  V"  _  (x  t  ..  (15) 


2 .  Second  order  body  correction 


The  velocity  induced  at 
excepting  the  :first  term  repi'esenting  the  source  j  itself  ,  is 


2  =  Xj  by  the  various  terms  of  (12), 


w, 


1 


fwi  .  !i  •  j  .  ,  f j_  ...l  +  fA .  i.  .  -f  i _ 

}  ik  2n  z-x.  2 n  2ilv  2n  x.-x,  2n  x.-x.-2ift 

■  J  zsx;  J  k  J  x 

■  e .  ie, 

+  -J-  m(-2ih)  +  W-&  ai(x.-x. -2ih) 

■  TI-  •  *■  „  IT  1  K' 


(16) 
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We  assume  in  the  sequel  that  x^.  is  an  increasing  sequence  and  that 

k  >  i,  With  £..  =  x,  -  x.  >  0  w.  becomes  . 
jk  k  j  J 

'  e.  .  fc.  .  £.  ,  ie.  ie, 

w.  =  ■r^-  ---.7-  -  \ —  ~ —  +  ij —  - — r^nrp  +  cu(-2ih)  +  - w(-liV-2ih) 

)  2n  2ih.  .  2n.  2n  £j^+2xh  it  v  n  jk 


By  the  sane  token,  the  velocity  induced  at  z  =  x^  is 


e.  -  e.  .  e .  .  ie,  ie . 

w.  =  ~  ™  ~  -  n1  - y,  —  ui(-2ih)  >  —3-  to(£  n:-2ih) 

k  2ir  2ih  2n  2ir  2-jjc'2ih  it  .  it  Jk 


To  satisfy  the  body  condition  (10)  we  have  to  cancel  and  v^  (16)  in 

the  vicinity  of  the  source  j  by  superimposing  on  the  original  source  an 
additional  source  of  strength  e.u./2n  =  0(c2)  and  a  vertical  doublet  of 
strength  -e4v^/2tr  =  0(e3).  Since  we  limit. the  present  computations  to  second 
order  effects,  we  may  disregard  at  this  stage  the  doublets  and  consider  the 
sources  solely. 

The  strength  of  the  additional  source  at  x.  is,  therefore,  given  by  (17) 

‘  •  -  '  :  ’  '  .  V 


"tor*.  4^  (-2,IO’2h)  +  4^  [- X~-  +  ~  ? H-l]k,-2h)]  (19) 


jk  >  - 


while  at  x,  through  (i 8)  and  (69)  we  have 

K_  -  • 


I--'3  4“7  (-2tro"2h)  -4;,T  It~  -  pT-HidT  +  2B(-£.t  -2h)-4ite  2hcos  (20). 


4iit  Lr.7'  ’  £?' +4h7 

jk  jk 


The  exprossionvof  the  stroamfunct ion,  associated  with  the  two  sources,  far 
downstream^, is  easily  found  from  (14),  (19)  and  (20)  as  follows 


ix.  ix 


\|/!*b(x,h)  = .  Im  e"ix{(c.^o  ^+e£c  +-ejcjJ(°  j-e  X) 


ixi  ixk  b 
>  J-0  K)il!  , 


:■ 


E'-. ' 


r 
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where  the  superscript  b  stands  for' body  corrections  and  the  coefficients  in 
(21)  have  the  following  expressions  , 


2e 


=  4e  ^cos  £.,  •  .  -  .  .  (24) 

•  •  jk  jk. 

b  - 

I:\j_  represents  the  interaction  between  the  local  terms  of  the  two  sources  and 
tends  to  zero  like  l/£..  ,  while  represents  the  interaction  between  the 

free  waves  trailing  behind  source  j  and  the  disturbance  of  source  k  . 

3,  Second  order  free-surface  correction 

The  free  waves  related  to  the  free-surface  correction  are  generated 
by  the  linearized  pressure  p**(x)  of  Iiq.  (6).  In  the  case  of  a  pair  of  sources 
p^  can  be  written  as  follows 

Pjk(x  *  .'jfrj’jJ  ir,tt  <irP)l  C  5 

the  different  terms  resulting  from  the  substitution  of  the  real  and  imaginary 
parts  of  wt.  and  wl.  (12)  into  (16).  Carrying  out  the  substitution  we 

JK  JK.Z 

obtain,  for  instance,  . 

PjkCx;xj -xk-h)  ■  ‘PVkVk'YkVdWWWVkVj1  (26) 

where,  for  brovity,  the  following  notations  have  been  used 


ai  =  a(x-x^,-h) 
2h(x-x^) 

*» 


3^  =  S(x-x£,-h) 


y£  3  '(x-x’)r+hY 


'£  (x-X^+h* 


(A*j,k) 
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is  the  pressure  resulting  from  the  nonlinear  interaction  between  the 
.velocity  components  induced  by  the  two  sources,  while  p^  (or  p^)  stems 
from  the  velocity  of  the  same  source  j  (or  k) .  It  is  readily  seen  that 


L  (x;x  h)  =  lim  |p 

•  JJ  J  x,  ->x .  J 

k  j 


(28) 


fkktexk'h)  * . yPjk 


x.->x, 
J  > 


It  turns  out,  therefore,  that  it  is  sufficient  to  determine  the  free  waves 
generated  by  the  pressure  distribution  p.,  and  to  take  the  limits  afterwards. 

JK 

The  complex  potential  representing  the  flow  generaled  by  p^  is 
(Wehausen  and  Laitone,  1960) 


fjk  »  -  “  /  Pjk(s;Xj,xk,h)  (D(z-s-ih)ds 


(29) 


If  p,,  ,  as  a  function  of  s  ,  is  absolutely  integrable  in  the  infinite  inter- 
val>  w  may  be  expanded  for  large  x  under  the  integral  sign  of  (29),  is 

not  integrable  as  it  stays  because  of  the  first  two  terms  of  (26).  We  are  going 
to  integrate  these  terms  first  to  eliminate  this  difficulty.  By  using  integration 
by  parts  and  (66) ,  (67)  wo  obtain 


Im  ~  /  [o.(sjx.,h)o.  (s;x.  ,h)  +  3 .(s;x .,h)e.  (s;x.  ,h)]w(x-s)ds  = 

it  j  j  K  k  j  J  x  K 


-i6ir2e'2hcosJljk  -  Im  ~  /  4i(B^Yk+YjBk+akXj+cxkX;j)w(x-s)ds 


(30) 


After  this  first  integration  tho  prossuro  is  intograblo  and  by  expanding  u(x-s) 
for  x  »  (Iiq.  65)  we  obtain  for  largo  x 
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#j,k  =  Im  fjk  -  -16n2e”  cosJLj.  +  Ira  8?,  /  [i (e^Y^+Yj ^k+ct j xH+ctkx j ^ 

•,i',kMA*Vj  *ii,k*V]*<)'i'*kBjlels,ls  <*->  (31) 

The  integration  of  the  various  terms  of  (31)  may  he  carried  out  by  using 
the  representations  of  a^,  6^  given  in  Appendix  (6,3).  We; have »  for  instance, 

_/"  (yk*ypolsds  =  fu.  f  Ts-^.lKH^pTiT  ll0,°iS<15  t32J 

I  l<yk*yk>«is<>»  .  i  /  [I!0  /  Jl,lclSlls  '3i» 

-00  *  *  .go  -  O  J  K 


and  by  using  also  (67) 


J  5j6kolsds  «  /  Yk  *j~  (^e^Jds  =  J  olb  7  dpUs  ' 


CP  ,  CP 


X  1  t  is . 

*  ff  i-W  as 


Using  the  residue  and  the  semi residue  theorems  and  the  definitions  of  a 

^  J' 

and  8  (63)  wo  arrive  after  integration  in  (34)  to  the  final  expression  of 

■*  .Is  _  i  y  iXl  ix.  ■■*,/■■• 

(x,hj  a  8 tie  Im  e  •  {(e  •'♦o  K)  (in/l  ♦  ~jr- ♦  a(-tu  0)  - 

Xjjk  JS, 

h  ix>  ix.  .1  2h 

-  2a(-£j^-2h)  +  +  (o  J-e  )ij-ttti^  'f~* 


vec-tjk(0)  *  2e(-tjk  -2nr  ♦  ♦ 

♦  2tfio  k{e  Jk+2e'21,e  *!*))  -  16a2o'?ht«)s£jk  (xH  (35) 
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i>.. and  iji  are  obtained  from  t|i..  by  taking  the  limit  i..  ■*  0 
J  3  KK  J  k  J  K 

(28)).  Carrying  out  the  computations  and  using  the  properties  of  the 
functions  or  and  8  (69),  yields  the  following  final  expression  for  the 

streamfunetion  far  downstream 


I  Is  -ix  ,  1Xi  ,  1Xk  s  s  1Xn  X\  s 

Jm  e  1X{ (eje  3+e*e  k)(As+iBs)  ♦  £jek[(e  J+e  k)C*. 


:"A'~  V ;  ijt .  ix  ix 

•-  r*  Ce;  ^ -  e  k)iE*k  +  e  ki(Ijk+iK?k)] }  +  (ej+e£)Cs  +  e^e^cosi^ 

*  '  •  ^  (x-*-)  (36) 


’.t  .X> 


where 


AS  *  /||^h[£n  2h  ♦  ~  ♦  C  -  2e'2h  Ei(2h)  ♦  ^  (37) 

VO  ~  ’ 

iSp  •  . 

2e’h(iA2o’2h)  (38) 


cs 

Ljk 

*•  ^  0:'k[4  tn(l+  |t~)  ♦  a(-£jk>0) 

-  2aC-tjki-2IO  ♦  **“^1 

(39) 

-  jk 

*  f  e%tan'1  ~~  ♦  B(-£,k  0)  + 

:f,.  .  v^,  £jk  .  Jk» 

2»(-.jkj-2h)  ♦ 

(40) 

Is 

iii. 

jj. 

(41) 

.,8/  i 
hk 

•<%/;-’4#%fi^or2h)|siii  .4^..  ’ 

(42) 

(43) 

The  superscript /  s  stands  for  f roe-surface  corrections,  C  »  0.5772  ...  in 
(37)  is  the  Kuler  constant  and  ifi(2h)  in  (37)  is  the  exponential  integral . 

The  first  term  in  (35),  associated  with  As  and  8s,  represents  the  free 

waves  generated  by  the  interaction  between  the  different  components  of  the 

velocity  of  thtr  some  source,  while  the  remaining  term,  in  c.e.,  is  related 

J. "  . 


HYDRONAUTICS,  Incorporated 


12 


to  the  interaction  between  the  two  sources.  The  last  term,  in  Gs,  is  the 
nonperiodic  D.C.  term  which  does  not  contribute  to  the  wave  resistance. 


4 •  Ajynjgtotic  expansions  for  the  free  waves  amplitude 

The  amplitude  and  the  phase  of  the  free  waves  depend  entirely  on 


£j»  ek»  Xj,  Xj.  and  the  coefficients  of  (36)-(42).  Since  wo  have  found  closed 
analytical  expressions  for  these  coefficients  it  is  easy  to  study  their  behavior 


for  large  and  small  Froude  numbers. 

Small  Froude  number .  We  consider  the  limit  h  =  h'g/U'2  •*■  «  for  a  fixed 
ratio  h/JLk  =  h * / Jl j k  ,  i.c.  we  lot  also  JLk  =  *jkg/U'2  «  . 

Under  thi..  imit  the  amplitude  of  the  first  order  waves  (15)  a*  is 
Q(eo  ,l).  The  amplitude  of  the  waves  resulting  from  the  bodv  correction  (22-24) 

j  rk  i. 

a 


IIl>  is  0(e2e  *Mjk)  or  0(e2o  3*1) ,  Finally,  the  largest  possible  amplitude 
of  the  waves  associated  with  the  free- surface  correction  alls  stems  from  As 


(37)  and  is  0(e2e’htn  h)>  Hence,  auu/al  »  0(n/a  )  or  0(eo-,il1)  and 


I  lb ,  I 


-2h, 


I  T  c  I 

a  /a  ■  0(etn  h). 


jkJ 


Since  e  is  the  inverse  of  the  thickness  Froude  number,  the  singular 


behavior  of  the  *>qo  waves  for  h  *  «  found  in  previous  works  (Salvesen  1969, 

Uagan  1971)  is  validated.  Moreover,  the  order  of  magnitude  of  a1 -/a1  for 

*  »  _  .  * 


1 1 S  I 

U-  *  0  is  exactly  found  a  /a  in  IM/W'2. 


This  point  is  discussed  in  more  detail  in  the  analysis  of  the  wave  resistance 
(Section  IV). 

h'n  »  .  ..  .  h 


Large  Froude  Kith  h  ■  |)“!f  k  O.  •»  o  and  --t  a  0(1)  we 

obtain  from  (15),  (23)  and  (39)  or  (40)  for  the  largest  possible  amplitudes 


I 


a  a  0(c), 


aIll>  « 


and  a*ls  »  0(«7/i -k)  . 


i|i 

Hence  the  ratio  a  7a  is  of  order  almost,  the  latter 


being  a  fixed  (generally  small)  number  for  a  given  body. 
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5 .  Free  waves  behind  an  arbitrary  source  distribution 

For  an  arbitrary  distribution  of  n  sources  of  strength  located 
at  x  =  x^,  y  =  0  (H=l,2, . . . ,n)  we  obtain  from  (15),  (21)  and  (36) 

I  -l  -i*  n  ix‘ 

i|»  (x,h)  =  Im[(-2ie  n)e  £  e.e  -1]  (44) 

j-1  J 


.II,  ..  .  lib  IIS  .  -ix,  V  2  rtixjr.s  b 

ip  (x,h)  »  i p  *  ip  =  Im  e  {  >  £<•  e  J  [A  +i(B  +B  )]  + 


j  =  l  J 


n-1  n 


ix.  ix 


ix.  ix, 


*  l  l  £,6.1(0  •’♦O  k)C*  ♦  (o  J-e  k)i(U?k*K?k)  ♦ 

j«l  ksj+1  J  K  JK  Jk  JK 

♦  elXk  i(I^I?k*iK,k)])  ♦  GS(  j>  ♦  2  "f  l  t  .t.cost  )  (45) 

JK  JV  j  =  l  J  .  j«i  k=j*l  3  K  Jk 


whore  tho  coefficients  of  (45)  are  given. in  (22-24)  and  (37-43). 


IV.  THE  WAVE  RESISTANCE 


1.  General  formulation 

The  wave  resistance  may  be  determined,  for  instance,  from  a  momentum 
balance  between  two  vertical  sections  far  ahead  and  behind  the  body. 

With  tho  free  waves  profile  (S,D)  given  by 

nl(x)  *»  -  *l(x.h)  -  -(v^eos  x  ♦  ^stu  x)  (46) 

n'Vx)  »  n‘2  “V*  *  *».' ^in8i,‘  x)  "  ist  (47) 


the  wave  resistance  at  second  order  (Salvesen, 


!-%!».) 


is  given  by 
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l)1  + 


I)11  ♦  0(e4) 


I2 

\h  + 

cos 


stn 


*  i  (♦.' 


II 


cosycos 


*  *  °<“)  <48> 


In  (46),  (47)  and  (48)  **QS,  *gin.  g»“g  and  ^[n  represent  the  periodical 

sin  x  .nd  cos  x  terras  of  (44)  and  (45) ,  respectively,  while  stands  for 

»  *  LOIiw  v 

the  last  term  of  (45).  ^  is  exactly  canceled  in  (47)  by  the  constant  part 

of  n1  such  that  the  average  free  surface  elevation  is  equal  to  h  for  x  •+  <*>. 
The  dimensionless  drag  in  (48)  is  defined  as  D  =  D'g/plJ'4.  Using  hqs.  (44) 
and  (45)  and  substituting  in  (48),  wo  obtain  after  some  manipulations 


,,I  -2h 

1)  “O 


n  n 


l  l  «« 

j  53  1  k=l  ^ 


ek  cos 


(49) 


,n  lib  .  i Is  _-h 
l)  =1)  *1)  .» 


n  n 


>b. 


0  ()!  Sinn..  ♦  (Bl+Ba)cosfc.k) 

)»l  k*>l  •  J  J 

m«l  j»l  k“j ♦ 1  J  ■ 


n  n-1  n 


‘  •  Kjk  *4».V)>- 


(50) 


fcqSk  (49)  and  (50)  give  in  a  closed  form  the  first  and  second  order  wave 
resistance,  respectively,  for  an  arbitrary  source  distribution.  and  i) 
the  nonlinear  body  and  free-surface  contributions,  may  he  computed  separately 
from  (50)  by  selecting  the  appropriate  coefficients. 

Kc  are  going  now  to  discuss  the  application  of  (49)  and  (50)  to  a  few 
simple,  but  instructive,  particular  cases. 
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2 .  The  isolated  source 

The  simplest  case  conceivable  is  that  of  an  isolated  source  beneath 
a  free-surface  (Fig.  2a)  representing,  at  second  order,  a  body  of  semi¬ 
infinite  length  with  a  parabolical  blunt  leading  edge. 

The  wave  resistance  components,  obtained  from  (49),  (50)  with  n=l, 
have  the  simple  expressions 


»l  = 

cV2h 

DIIb* 

e3(-2o'4h) 

(51) 

uIIS  « 

eJ[-2e  (l+2e  )] 

The  wave  resistance  of  (51)  is  given  in  Fig.  2c  by  using,  however,  the  more 
common  drag  coefficient 


Hence,  with 


we  have 


•Cjj  •  ».  nyq.W*r.  and  the  parameter  T'/li'  »  t. 


»  t(‘  ^  4  t2!’1  ^ 

■V  '  li  5  4*|» 


C*  "•  3lte‘a‘ 


c‘,b  •  -4IA"’11 


(S?) 


c;,,s .  -4»-ai  *<..*•*) 


Ute  first  order  coefficient  C*  as  function  of  li*//gh'  =  i/VK  has  the  well- 
known  shape  with  a  maximum  at  U'8/gh'  »  2.  The  second  order  resistance 
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coefficients  are  negative,  i.e.  they  always  diminish  the  values  cf  the  first 
order  drag.  To  illustrate  the  relative  magnitude  of  the  second  order  effect, 
the  ratio  /('  ,  which  does  not  depend  on  t  ,  lias  also  been  represented 

in  Fig.  2.  For  U'//gb'  <*  this  ratio  tends  to  zero  like  -8h  and, 

therefore  -*■  -8ht  =  -8c,  Hence,  1/e  =  U'2/gi"  has  to  be  much  larger 

than  8  in  order  to  ensure  small  nonlinear  effect  at  high  Froude  numbers.  For 
U'/i/gh'  0,  »  -2h  i.e.  t2C|,)1/t('jI)  -*•  -2e,  and  the  nonlinear  freo- 

surface  correction  becomes  much  larger  than  the  first  order  wave  resistance 
because  of  the  f roe-surface  effect.  In  the  intermediate  range  C^/C^  grows 
monotonicully  (Fig,  2). 


Finally,  to  estimate  the  relative  importance  of  the  two  nonlinear  effects, 
the  ratio  !  |)lls/n^  s  is  also  represented  in  Fig.  2c. 

For  *1 ' *  ■*  "  this  ratio  tends  to  3.  As  li'//gh*  decreases  1»  /U,  tends 

to  infinity  like  e“  ,  lienee,  in  the  present  ease,  the  tree-surface  correction 
is  always  much  larger  than  the  body  c  rrect ion . 


The  influence  of  the  flnctvtas  of  the  leading  edge  shape 

To  study  the  influence  of  the  shape  of  the  leading  edge  we  have 
cons  idcrcU  a  semi  - 1  «f  in « to 'body  generated  by  ten  sources  of  e«pta  1  st  rength 
and  at  equal-  spacing*,  i.e.  w  *  t.  9  t/lfl  sy  L  * .’(k»j)h/D  (Fig.  21*). 

.  j  *■>  )  k 

This  distribution  represents  quite  accurately  a  wedge  shaped  nose.  It 
turns  out,  now,  that  the  general  drag  fomtlae  (4D)  and  (at))  simplify  in  this 
case  and  become 


Si*  *  ltf~«V~lj 


10  to 


j*i  Mi 


cost  , 

JW 


m 


ouu  - 


.  ,,  lo  10  10  *5  10  . 

•V .  j  -h,J»  v  V-  -  t  r  c  ,o 


■t  .1  cos*it  ♦. ! .  I  j  _  iu«sy  * 


jif  h«l 

10  It) 


k5t  ,  l  .  I-  J.  ’:V' 

}  *3*1  j*U*j*l 


*  . 1  .ah  »"  **’  to  o  lo 

**  cse'  ‘(2c“  )|  |  j>  cost  ,  *2  J  f  cost  .cost  | 

j*i  Mi  ,tt*i  j*i  i*}*i  "*  }h 


<S4) 
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I  Is  -l.  s  10  10  10  9  1° 

U  =  -10  t3e  cosijk  *  l  l  l  (i*{CQS*  .-K^.sint  .  )]  = 

j  =  l  k=l  JK  ml  j=l  k=j*l  ->k  mk  ->k  mK 


•10~’>t3[2e  ^lCl*2c'2*')  f  cost..  ♦  2o*2h(l+2e'2h)  • 

j=l  k=l  JK 


10  9  10  .  10  0  10 

l  l  l  cost  cos^  -  4e'-h(l-2e'2h)  £  \  f  sin*  .  sint.k} 

m=l  j=i  k*j*l  wk  J  c*l  j=l  k*j+l  aK  Jk 


0*/e2,  l)*VVc-3-  and  t)us/e3  have  been  computed  as  function  of  h  for 
different  ratios  between  the  fore  body  length  and  the  suhaergence  depth  (h*/h') 

In  big.  2c  we  give  the  results,  again  in  terms  of  t‘  *  nVfl.SoU'^T  « 

I  -I[  w 

■  tC,,  ♦  t?C(>  and  t  «  r/h*.  for  hW  -  9  ami  k’/h*  «  4. a.  The  first 

order  resistance  coefficient  C*  drops  markedly  as  the  ratio  h‘/V  increases. 

lor  liV*/gh‘  *  3  the  curves  beeeae  very  close,  because  the  ratio  between  the 

length  of  the  free  waves  x*  *  2*M#tyg  and  the  length  oi  the  fore  hedy  V  is 

larger  than  2*  for  the  largest  h*  considered,  At  smaller  t  ruudc  mcahcrs, 

however*  the  interference  effects  are  esahifest, 

A  sore  interesting  point  Is  that  related  to  the  influence  of  the  fineness 
on  the  ratio  ,  (l  Ig.  U*  >*eepiihg  for.  iarge..  (U*/#iT  *  3)  orMll 

(U'/^gh*  *  0*7)  broude  mashers,  the  relative  nonlinear  .effect' is  smaller  for 
fine  shapes,  as  much  as  half  of  that  of  the  blunt  shape  (for  j.'/li'  *  i»).  for 
the  latter  case  is. almost  constant,  for  3>:>  tiygh*  >  1.7  and  approx i- 

stately  equal  to  »V?$.  . 

•  finally,  the fifignesa  has;  little,  influence  on  the  ratio  i>* *s/!}* 

(fig,  2c)  which  reouius  large,  irrespective  of  the  shape  of  the  leading  edge. 


,Iis 
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4.  I nt ert'ercnce  cf f cct s ;  the  source-sink  body 

We  consider  now  a  body  generated  by  a  source  and  a  sink  (big.  3a), 
which  is  the  simplest  case  of  a  closed  body. 

Using  again  the  general  expression  of  the  drag  ((40,(50))  with  =  c, 
.Cj  =  *£  and  tp  55  1.  we  obtain  in  this  case 

l>*  »  e^2e  “^(1-cos  I.) 


,IU>  i  -hr  h  ..  .  . .  .  b 


>c3e  ’(-2hp(l-cos  i.)  ♦  Ip (1 -cos  h)) 


(57) 


Ulls  « >eV4[2Assin  h-dt  j^sin  l^M.|,(l»cos  !.)♦! p£l-eos  W)*kpSin  I.)  (SS) 


the  different  coefficients  Wing  given  by  (22)- £24}  and  (37) - (42) .  We  have 
selected  an  example  with,  h*  «  2oh*  consistent  with  considering  elongated 
bodies  reselling  ships.  In  this  c«.  v  we  ecu  use  the  asymptotic  expressions 


of  a  and  8  which  appear  in  ft,  and  1  * ,  (30)  ,(4i>g  to  obtain 

*.  .  ’  ■  In  '  ■  .  1  H 


;a- 


2  ;  *h  !~Sh  •  2h* 


-1  3-  2h  .  .-1  ' " 

lij^,  *  ••  U(t  )  .  • 

these 'coefficients- being  related  to  iW  interaction  bet ween. the  local  dlstdr* 

a  ,  | 

Ua»«e»«  Ihtiy  terns  of  order •  l  •»  (2«hV  or  larger  have  been  retained  in 
($$)•  Id  plot  the  different  drag  .eo^hment ,  we  have  used  thoiwue  -faint iar 
di«0«si©niess  jwraifteter^  Cj  ®  U*/«d3*sh*  -  and  e(  *' 1' /!,*  .'  -With  this  r«pre- 
sentat ion  we  have  - 


h  H  h.  ‘i.  \ 
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where 


T  t  t  r 

~  hi)  =  2()hl) 


(p)2  h2l) 1 1  =  400h2Dn 


fin.  Fig..,  3b :  represented  separately  as  functions  of 

th^ev;.^)Qct^<^s  -ha^hg  an  oscillatory  character.,  hike  in  the 
previous  cases,  the  body  correction  is  smaller  than  the  -freo-surface  second 
order  correction.  For  sufficiently  iiigli  Rroude  numbers  (U'/'^b1  >  0,3) 
is  roughly  0.2  »  0.3C**,  while;  for :  3  _  its  relative 

magnitude  becomes  oven  smaller,  Moreover,  for  the  ratio  '  L'/h'.  =  20  considered 
haro/the  local  verms  in  (5?)  and  (38),  i.r,  the  terms  in  l:^,  As,  C^,  K®^ 
are  small  compared  with  the  other  terms  and;  for  the  sake  of  the  analysis, 
and  1)^®.-  may.  be .  approximated  for  ,0.15  ;  U 1  /  /gLV; 1  d .  45  ...  by 


-  e3o’*‘l^cos  1.(1 -cos  I.)  a  -  e^c'^cos  L(l-cos  L) 


l)1  -  *-6V2hl®.,(WcoS  I.)  634e'2h(l+2e‘2h)cos  1(1 -cos  L) 


/  Hence, b*  (Sf)  .oscillates  like  1-cos  1  while  l)^s  and  0*^.  bohavo 
"'like  cos:  1{1  cos  b).  As  a  result, the  second  order  correction  tends  to 
sharpen 'the  peak's  oi  the  resistance  curve  and  to  make  the  hollows  wider.  At 
'the  .peak. of  fl  2  resistance  curve  for  ll'//gb'  *  0.32  the  ratio  I)11/!'1  is 
i  of  order  40 and  this  ratio  grows  as  the  Froudo  number  decreases, 
like  in  the  case  of  the  semi-infinite  body,  the  freo-surface  second  order 
effect  becomes  unboundedly  larger  than  the  first  order  drag  for  U’/Zgl.’  •+  0, 
dominant  tern  In  {58)  being  A  *v  o  nn  h. 
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V.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS  ;? 


Although  an  useful  closed  solution  for  the  nonlinear  wave  resistance  in 
two-dimensions  has  been  obtained,  the  main  purpose  of  the  present  work  is  to/ 
draw  qualitative  conclusions  on  nonlinear  ship  wave  resistance.  Two  features. 


of  the  results,  summarized  in  Figs.  2  and  -3,  seem  to  be  general:  '(L)  the  free- 


surface  second  order  contribution  to  the  drag  is  larger  than  that,  of/the.  body 
correction,  and  (ii)  at  small  Froude  numbers  the  nonlinear  effects becqrae/yery; 
large  and  the  basic  expansion  of  the  velocity  field  is  singular. 


The  results  depicted  in  Fig.  2  seem  at  first  glance  of  limited  interest,'  v 
because  theyfappiy  to  open  bodies .  But  there  are  indications  that  due  to  • 
viscous? effects  the  stem  waves  of  ships  are  smaller  than  those  predicted  by 
the  sink  distribution  representing  the  geometry  in  the  linearized  theory.  Such 
an  assumption  is  tantamount  to  considering  an  open  body  so  far  as  wave  resistance 
is  concerned,. /The" results  of  Fig.  ,2  may,  therefore,  have  some  bearing  on  the 
part  of  the  wave  resistance  Associated  with  the  bow  waves.  It  suggests  that:  . 

(IX  the  •tionlinedr  wavef resistance  is  sinutfer  than  the  first  order  drag,  and 


'■J '  ‘  \  - 


(ii.)  that  the . second/ordor  effect  is  less:  important  for  fine  shapes  "than  for 


full  fornix 


I.'  is  the 


s  of  rth.e  't'of  e  bouy.  L‘or;  !;ipdttrate  i:roude  nuribcys-/  (U  P  ‘A ,  ..where 

e  fore  body  length)  .  ,  -//'  . ; 


The  .interference  effects  botwconther  lead  inland  trailing  edges  /are;"-:, 
reflected  by  the  results  of  fig.  5,  Obviously, .  these:  effects  -Will 'he  largely. 
reduced  if  the  trailing /edge, is  assumed /t  a.  ;be -.of  a/fino  .shape  or  open,  (to;/:/- 
account  for  wake,  for  instance} ,.  At  any  rate,  tlie  present  computation.S“indi- 
catc  that  tho-'re'sistanco-pcahs^arv/^.jp'ion-tedXiy'.rsccqnd  order  effects  while  the  . 
hollows  become -wider.  ••  :  •/•/•" 


The  method  used  in  tho  presciit  work  may. /be  easily  extondeA'ft.o' account  Jpr 
vorticity  in  two  diwepsions,'  It  ogtabiishet  a  pattern  which  "maybe  useful;! n  . 
attaching,  the  problem  of  throe-d iwons ionai  ^eiwrcc/^ncrato4.J>pd-lca'. 
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VI.  APPENDIX  -  The  Properties  of  the  Function  u>(s) . 


The  function  w(0  =  ct(£,n)  +  i$U>n)  is  defined  by  (13)  as 

V  '  -  >  01(0  =  .  e“li:Ei(iO  .  =  -e'1C  /C  £dX 


in- the  X  complex  plane  (Pig.  4)..  Another  expression. 

obtained  by  the  change 

of  variables  p  =»  -X +  j,  is 

.  '  ■■  V- 

y_  '  -  ; 

'  '  e"ip 
”u)  ■ 

(63) 

with  the  p  plane  depicted  in  Pig.  4. 

oi  has  the  following  asymptotic  expansions: 

ir  “  m 

01  a  e“  •  (iit  +  C  +  tn  C  +  l  )  • 

m-1  “■"1  ■ 
and 

(t  +  0) 

(64) 

v  (n-l)i  , 

u  *  l  n  +  uJO 

m=l  (i;) 

(C  *  *) 

(65) 

where  oi  (C)  *  2trie"15  for  Re  q  +  »  and  u  »•  0  for  Re  i;  *  . 

00  00 

series  of  (65)  is  divergent. 

The 

Herewith  a  few  relationships  used  in  the  present  work: 


(i) 

and  particularly 


do) 


uU.n) 


6U,n) 


-  iw  +  — 
* 


36  n 

3?  ■  }?~~+~nr 


Du  i 

or.  ‘  F"V”nT 


(06) 


(67) 
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2Trie-1^  +  +  in) 

a(ij,n)  =  2ireri  sin£  +  a(-5,n) 
B(S,n)  -  2Tten  cos£  -  B(-C.n)  .  ■ 


(ii)  w(0  = 

which  leads  to 


(68) 


(69) 


To  compute  a  and  8  one  may  use  the  series  of  (64)  for  |s|  <10 
or  series  (65)  for  |sj  >  10.  In  the  latter  case  the  series  has  to  be  truncated 
at  m  =  N  where  N  ^  |e|  . 

u(e)  is  related  to  the  function  ii, (?)  tabulated  by  Abramowitz  and 
Stegun  (1964)  as  follows: 


wU)  «*.  e“1C  [2iti  -  E^-i”'.] 


(70) 
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•INFINITE  LENGTH:  (o)  A  SOURCE  GENERATED  BODY, 
i  ( c  }  WAVE  RESISTANCE  CURVES . 
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FIGURE  4  -  THE  A  AND  p  PLANES 


